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The renewable energy sources (RESs) projects are solutions with 
environmental benefits that are changing the traditional power system 
operation and concept. Transient stability analysis has opened new research 
trends to guarantee a secure operation high penetration. Problems such as 
frequency fluctuations, decoupling between generator angular speed, 
network frequency fluctuation and kinetic energy storing absence are the 
main non-conventional RESs penetration in power systems. This paper 
analyzes short-circuit influence on frequency response, focusing on weak 
distribution networks and isolated, to demonstrate relevance in frequency 
stability. A study case considered a generation outage and a load input to 
analyze frequency response. The paper compares frequency response during 
a generation outage with a short-circuit occurrence. In addition, modular 
value and angle generator terminal voltage affectation by electric arc and 
network ratio R’X, failure type influence in power delivered behavior, 
considering fault location, arc resistance and load. The arc resistance is 
defined as an added resistance that appears during failure and influences 
voltage modulus and angle value results showing that intermittent non- 
conventional RES participation can lead to frequency fluctuations. Results 
showed that arc resistance, type of failure, location and loadability determine 
the influence of frequency response factors in weak power systems. 
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1. INTRODUCTION 


Since the last decade, renewable energy sources (RESs) penetration is changing the electric power 
systems (EPSs) energy matrix. Agents used to reduce conventional unit participation in daily load graph 
coverage [1]-[4]. RESs bring multiple benefits for EPSs and environment, but some of them, such as 
photovoltaic and wind power generators, have intermittent electrical energy production, which can produce 


frequency fluctuations [5]-[11]. 


Another disadvantage is that RESs connection to the EPSs is through converters, this produces 
decoupling between generator angular speed and network frequency, as occurs in several wind turbine 
applications [12]-[18]. Also, solar photovoltaic generation is static and do not store kinetic energy. These 
two technologies do not contribute to primary frequency control [19], [20]. Static converters photovoltaic, 
wind generation and batteries energy storage systems affect the total EPSs inertia ratio, so this effect on 


Journal homepage: http://ijece.iaescore.com 


1078 O ISSN: 2088-8708 


frequency deviations has become a topic of interest when an EPSs consider the participation of these 
electrical energy sources with high penetration [21]. Frequency stability emerges as the major problem in 
modern EPSs operation. The frequency variation in the EPSs voltage is one of the most important operation 
variables. Its magnitude deviation, regarding its nominal value, can have detrimental effects for both 
consumers and EPSs elements that make them up. Frequency deviations are inversely proportional to rotating 
total inertia of the machine and system disturbance caused by them. This disturbance can be a load variation, 
generating unit disconnection, transformer or transmission line opening, or a short circuit event. 

In this context, isolated EPSs show a more pronounced event behavior of this phenomenon since 
they are weak systems that do not have interconnection with other areas allowing frequency behavior 
improvement against disturbances [22]—[24]. Also, they have a lower inertia and supported by RESs. Many 
researches in this field aim to address this problem and its solution. Unlike isolated systems, large EPSs has a 
large-scale generation units’ arrangement, they are robust interconnected transmission networks and operate 
with spinning reserve values which guarantee that frequency stability is not a dangerous phenomenon. In 
these EPSs, the attention to voltage collapse and its related problems is a priority [25]—[27], transient angle 
instability problems and a recent priority search for characterize small disturbance angle instability problems. 

CIGRE/IEEE defines different EPSs stability categories [28]. This definition shows each category 
and analyzes the procedures considered on its evaluation. RESs in isolated EPSs evaluating methodologies 
are detailed in [29]. Especially in frequency stability studies, emphasis is placed on failures to be analyzed as 
well as on generation units’ behavior and EPSs in general [30]. In study [29], it is described the behavior of a 
generator during faults. They speed up, followed by a sudden load reduction in voltage collapse caused by 
faults; the rotor angle also experiences an increase when the machine acceleration occurs. This behavior is in 
correspondence with a frequency increase. While a short circuit effect on a high nominal voltage power 
system occurs, control agents used to notice a frequency increase most time caused by the failure event, then 
a voltage drop is reported and there is a resulting link between machines, showing a high X/R ratio. 

However, generalizing this is incorrect, because of there is a dependence on this behavior with 
machines’ load initial state and X/R relationships that may remain in network equivalent positive sequence 
links. Isolated local networks have specific characteristics and favorable conditions can be present in order to 
failures cause a frequency reduction that leads low-frequency and/or frequency variation protection 
operation. Real cases in which ground faults have caused a frequency decrease and low frequency 
disconnections have presented. This has been this research motivation which aims to analyze short-circuit 
influence on frequency behavior in isolated electrical networks and to show this phenomenon analysis 
relevance in frequency stability. A sudden frequency reduction in short circuit events and low frequency 
protections operation make necessary to include failure simulation in frequency stability studies. These 
aspects impact the generation penetration effects with static converters. 


2. RESEARCH METHOD 
2.1. Frequency stability in electric power systems (EPSs) 

Modern literature defines the frequency stability in an EPSs, as the ability to maintain an adequate 
frequency value in a steady state caused by sudden generation-load imbalance occurrences [28]. Instability 
events include sustained frequency oscillation that can cause generators and/or loads to trip. Four well- 
defined stages appear in frequency behavior caused by a generation and load imbalance. Figure | show these 
stages. The first stage is the inertial response of the units, in which the kinetic energy stored in the rotating 
units assumes part of power imbalance and therefore contains frequency deviation. Second stage is the 
primary frequency control, in which the speed drive of the synchronous units begins with the aim of 
frequency restoring and stabilize into a value close to nominal. The frequency steady state error is the 
difference between nominal frequency value and the value at which frequency reaches its stable levels, this is 
called frequency steady state error. The minimum frequency value in EPS response is conditioned by these 
two previous stages. This point is relevant to EPS frequency stability. 

Subsequently, a secondary control action occurs, which allows frequency recuperation to its nominal 
value minutes after disturbance occurrence. The fourth stage is the tertiary control, which depends on EPSs 
generation matrix, as it depends on the possibility to have a quick-start unit that allows the restoration of the 
network secondary control reserves. A unit generation loss or a sudden load increase cause frequency 
behavior described in Figure 1. During short-circuits generators usually speed up, causing an increase in 
frequency equations that describe synchronous generator behavior. 


2.2. Synchronous generator behavior 


Two opposite torques causes the synchronous generator rotation. The mechanical contribution 
caused by turbine Tmand the electromagnetic torque Tę. Mechanical torque speed up the rotor, while 
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electromagnetic torque slows it down. In normal operation these two torques are equal, therefore there is no 
acceleration torque T, and the synchronous generator rotates at synchronous speed. Equation (1) describes 
this behavior: 


dor 
I ype = Ta = mn- Te) (1) 


being J the inertia moment and 6, the rotor angle. 
Typically, this rotor motion is expressed in power terms instead of torque, and as an angle 6 
function that rotates at synchronous speed, therefore as shown in (1) can be written: 


2H d*6 
w, dt? (2) 
where: w,: synchronous speed 
H : Inertia constant 
Pa: Acceleration power 
P n: Mechanical power supplied by primary source. 
P,: Electrical output power 

Equation (2) defines the synchronous machine oscillation equation. An important term in EPSs 
stability, H constant or inertia constant which defines the machine capacity or kinetic energy and as shown in 
(3) defines it as follows: 


H= Exinetic _ 1 Jws 
Smachine 2 Smachine (3) 


where Smachine represents MVA machine capacity. This constant appears within machines defined data to 
carry out stability studies. The stored kinetic energy provides inertial response and primary frequency 
control. 

To determine the electrical behavior of the synchronous generator (4)-(9) calculate the machine fem 
variation in direct axes and in excitation field quadrature, all expressed in per unit on nominal data of 
machine capacity, frequency, and voltage. 


Frequency 
A 


Time 


Inertial Response Secondary Frequency 


| 
| 
| Primary Frequency Control 
| Control 
bosada oo 
frn=Nodir 
——— pe y 
Seconds Minutes 


Figure 1. Frequency behavior stages in generation loss event 
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dE; 1 i : 

—i= = eg = (Xa = Xaia E E] 

dt Ta (8) 
dE, 1 , 

—— = — [E} — (XL — X,)i, — Ea] 

dt To d ( q a) q d (9) 


In these equations E4, E4, E; and E% are sub-transient and transient fem corresponding to direct axes 
and machine field quadrature; V4, V,, ig, and i, terminal voltage and armature current components in both 
axes; Ey armature fem proportional to machine applied excitation voltage, Tios Taos Tao and Tj both axes 
sub-transient and transient time constants corresponding to buffering and field windings; Xj, X}, X 4 and X4 
quadrature and direct axis sub transient and transient reactance and Rą machine armature resistance. That 
variables and sub-transient constants establish a relation with buffering and field transient winding flux 
concatenations. In these equations when applied to protruding pole machines, if there is no buffering winding 
in interpolate area or its effect is neglected, a necessary adjustment must be made in corresponding direct axis 
equations since E4 does not exist under these conditions. The (10) and (11) calculate active and reactive 
power delivered by machine at its terminals, respectively. 


Puetiverea = Vatg + Vala (10) 
Qaetivered = Vzla ~ Vaig (11) 


Magnitudes transformation on d,q axes to corresponding magnitudes on real and imaginary axes 
can be calculated using (12) and (13) relation: 


kino cosa a 
H = ea sin(d) | F] 
Va sin(d) —cos(d)} LV; (13) 


where 6 is rotor angle, as expressed in (2), measured with respect to a real axis that rotates at synchronous 
speed and that serves as a reference to all system phasors in which machine is connected. This angle is 
calculated for each machine using (14). 


|Egle’® = Vin + (Ra + iXq)lrm (14) 


Equation set from (1) to (9) allows to obtain machine rotor speed and position while known its terminal 
voltage, which must be calculated from system electrical network solution by using load flow equations. For 
this, nonlinear differential and algebraic equations solving numerical methods can be used, which allow 
calculating and performing system voltage and frequency behavior in disturbances presence (15) is obtained 
from (4) and (5) which means currents can be calculated as follows: 


ilare elev 
ial R24+X)Xq E,—Va (15) 


Substituting in (10), power delivered is expressed in (16) as a sub-transient fem and terminal voltage 
function: 


Paetivered = ( Ze EV, cos(d = Oy -0 ad 


RZ + XX% eee 
+ ZE Vm sin(ô — Oy — 04) — Rav, + 5 Xa 
— Xq)Vnsin(2(6 — Oy))] (16) 


This (16), as a variable and sub-transient constant’s function, can be expressed in the same way with transient 
and steady state variables. In this case are required to adjust according to machine field type and consider that 
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in stable state fem direct axis is zero. Comparison between one machine power delivered and those obtained 
by a general numerical method in presence of system disturbance is permitted by this equation. In 
asymmetric fault case, the following correction must be made by using (17). 


Pactivedetivered = Puetiverea — [R21 + (Ro + 3Ry IG] (17) 


Power (P,) to be used in (2) will be given by P4eliverea OT Pactive delivered according to each case, plus 
armature copper losses, given by Relig + iz) = R,I?. Where Å, I, and lọ are machine phase currents 
symmetrical components. 

All the (1)-(17) can be used in iterative calculation processes if saturation effect is considered. The 
(16) is exact for initial pre-fault and fault conditions at the time instant in which failure occurs. This equation 
graph as ô function could also be obtained by an iterative process to readjust saturation values. 


2.3. Isolated electrical power systems frequency behavior 

Isolated electrical systems have characteristics that make frequency stability an important aspect. 
These electrical systems power plants are usually supported by diesel internal combustion engines. These 
units have the particularity of low H constant. That connected diesel units cover high demand value is 
another aspect that characterizes isolated EPSs and therefore when a generation trip occurs, an appreciable 
load fraction of total generation is lost. Consequently, isolated electrical power systems when face generation 
loss or sudden load connection, frequency fall rate is high and often low-frequency protections operation 


and/or df / dt is necessary to avoid EPS collapse. This is accentuated with RESs insertion with technologies 


decoupled from network through converters, which makes EPS further reduce energy stored in synchronous 
units. 


3. RESULTS AND DISCUSSION 
3.1. Study case electrical power system description 

The selected study case is an isolated EPS which has an internal combustion units’ generation 
composition, with low inertia values. During daylight the attendant units are photovoltaic power plants, 
located in two busbar sand a wind farm in five. There are seven loads located, as is shown in Figure 2. 
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Figure 2. Electrical power system studied 


Residential load predominates in the system, with 17 MW of maximum demand approximately 
during the first night hours. Figure 3(a) shows electricity demand behavior when the maximum demand 
reached in 2019. Figure 3(b) shows the typical behavior for a solar and wind generation summer day. 
Figures 3(a) and 3(b) analysis exhibit the moments where demand peak is highest, which corresponds to solar 
photovoltaic generation reduction with and wind generation poor contribution. This is an interest moment for 
frequency regulation. 
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Figure 3. Demand and RESs generation profile where (a) demand daily behavior and (b) typical solar and 
wind generation in electrical system study case 


3.2. Real event modeling 

One of the most important aspects in electrical power systems studies is to establish an adequate 
correspondence between simulations and actual EPSs behavior. Results precision, conclusion, 
recommendations, and measures proposed pertinence in carrying studies out to guarantee previous aspect. 
Following modeling correspond to biphasic ground fault occurrence in an electrical system line of the study 
case. Unsuccessfully protections reclosing occurs as failure is permanent. Second protection operation 
performs permanent opening in the failed line. 

Network operation state corresponds to 12.5 MW of demand. Table 1 shows voltage, active and 
reactive power demand on each load Vo, Pp and Qy, respectively. In generation, nodes case voltage and 
active and reactive power generation Vp, P and Qg are shown. 


Table 1. Electric power system initial operating conditions 


Bus Vo (pu) Po (MW) Qo (Mvar) Po (MW) Qc (Mvar) 
1(Two Units) 1.04 5.4 2.04 6.26 4.32 
2(Two Units) 1.00 0.0 0.0 6.26 2.34 

3 0.98 1.49 0.78 0.0 0.0 
4 0.98 0.0 0.0 0.0 0.0 
5 0.97 0.0 0.0 0.0 0.0 
6 0.97 1.02 0.70 0.05 0.0 
7 0.98 2.37 1.23 0.0 0.0 
8 0.96 0.57 0.19 0.0 0.0 
9 0.96 1.25 0.70 0.0 0.0 
10 0.96 0.40 0.40 0.0 0.0 
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Figure 4 describes the frequency behavior after disturbance occurred. The behavior is in 
correspondence to frequency meter description that recorded the event. Note that in moments after failure, 
frequency drops until protections operate first time, frequency recovers. First reclosing moment occurs, and 
failure event remains, a sudden frequency decrease occurs, causing low-frequency protections operation. 
Subsequently, definitive protections operation occurs. 

This event was modelled using data obtained from utilities generators, speed governors and voltage 
regulators modeling were carried out in correspondence with manufacturer data offered. Figure 5 shows 
frequency behavior, the minimum frequency value obtained as well as disconnected load because of low 
frequency are in correspondence with what happened in the actual event. The Accuracy level reached help to 
establish an adequate correspondence between modeling and actual system. 


62 
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00:34 00:35 00:36 00:37 00:38 


Figure 4. Actual frequency measured by frequency counter during failure studied 
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Figure 5. Frequency vs time in failure simulation 


3.3. Generation loss vs short circuit 

Frequency stability studies usually test frequency behavior in generation loss event or sudden load 
input. As shown in actual fault analyzed, short circuits effect in this network type can cause frequency 
instability problems. A comparison between EPSs frequency behavior generation loss event and short circuit 
occurrence is shown in Figure 6. During simulations system demand was 15 MW and a 3 MW of generation 
was disconnected, which represents 20% of total system demand, and the short circuit analyzed was a single- 
phase to ground at node 3 with a 19 Q arc resistance. 

Figure 6 describes the frequency behavior when these two faults occur. Note that in short circuit 
occurrence case, minimum frequency value (nadir) is less by almost 4 Hz than the case of 20% generation 
loss. In both cases, low frequency protections operate, removing 5.8 MW of load from service. 
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Figure 6. Frequency behavior in short circuit and generation loss events 


3.4. Generation behavior affecting factors 

That independent variables of generator constructive characteristics influence its behavior in short 
circuits events is allowed by (6). Generation power supply is influenced by the network place where failure 
occurs. The presence of electric arc and network ratio R/X are factors that affect modular value and angle 
generator terminal voltage. It also influences electrical power and internal generator fem in pre-fault 


conditions. Failure type is another factor that influences in generator power delivered behavior, as can be 
seen in (17). 


3.5. Fault location 

The failure location influences network frequency behavior. A comparison between frequency 
behavior when a single-phase short circuit to ground occurs at node 2 and at node 9 is shown in Figure 7(a). 
In both cases, when short-circuit occurs, frequency falls unlike what is traditionally expressed stability books 
of electric power system. In large electrical power systems, power supply by generators decreases however, 
in both cases there is an increase in electrical power delivered by generator 1 during failure as is showed in 
Figure 7(b). Frequency falls when a failure occurs at fault at node 2 and 9 but there are substantial differences 
in the minimum frequency value, and in frequency value rate change. 
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Figure 7. Short-circuit analysis (a) frequency behavior over time and (b) electric power delivered 


3.6. Failure type 
The failure is another aspect that influences generator behavior. Figures 8 shows the frequency 
fluctuations for different short circuits type at node. That failure occurs with a 19 Q resistance electric arc, 
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and 150 ms failure duration is considered for several cases. There are differences in frequency behavior prior 
to different types of short circuits occurrence for single phase to ground and biphasic to ground frequency 
drops. In a single-phase short circuit to ground case, a low-frequency protections operation occurs 420 ms 
after the occurrence of the failure, removing 5.8 MW from service. In biphasic short circuit to ground case, a 
minimum frequency value does not reach low frequency protections settings and therefore, in failure event, 
there is no effect in users. 

Separate mention deserves three-phase short circuit to ground. Initially, frequency increases until 
protections act and clear the fault. At this moment there is a sudden frequency decrease to values below low 
frequency protections adjustment, causing a disconnection of 5.8 MW, 1.56 seconds after fault occurrence. 
Either in single-phase short-circuit to ground case or in three-phase to ground case that frequency instability 
phenomenon occurs Figure 8. 
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Figure 8. Frequency behavior for different types of short circuits 


3.7. Arc resistance influence 

The aspects tested above the impedance influence by generators when the fault occurs, and the 
terminal voltage imposed by the network on the generation nodes. The arc resistance is an additional 
resistance that appears during the fault and influences the value of the voltage module and angle. Figure 9 
shows the behavior of the frequency when a single-phase short circuit with arc resistance occurs and when it 
is clean. 

Figure 9 shows the difference in the frequency behavior when for the same type of short circuit and 
in the same node this occurs through arc resistance and when it is cleared a short circuit. The minimum 
frequency value in the case study occurs with an arc falls 1.21 Hz more than in the case where the fault is 
cleared. In both cases, the frequency drops, however, in the case where it occurs through an arc, the low- 
frequency protections act, removing 5.8 MW from service. Therefore, it can be considered that there is a 
frequency instability problem in the case of the arcing fault, and yet the non-arcing fault itself does not cause 
an instability problem. 


3.8. Effect of the renewable energy penetration 

The frequency behavior in electrical systems is directly related toto the kinetic energy stored value 
in synchronous units that is certainly defined as real inertia. The photovoltaic penetration in the SEPs 
replaces synchronous generation. Also, for the same demand value, there is less inertia as the solar 
photovoltaic power increases online. Figure 10 shows the frequency behavior for two different generation 
compositions when a biphasic short-circuit to ground occurs in the power line that interconnects nodes 9 and 
10, with a cleaning time of 350 ms. For a demand of 12.5 MW, two cases were analyzed, the first without 
photovoltaic solar generation and the second substituting a generation unit for 2.8 MW supported by solar 
generation on the grid which represents a 22.5% penetration. 

Figure 10 describes a marked difference in the frequency drops in both cases, the inclusion of power 
units that do not present inertia causes a 2.5 Hz additional frequency drop, the low frequency relay is 
activated avoiding an instability problem. Also, it is noticed in this simulation that with a higher renewable 
energy penetration exists a great speed frequency variation because the inertial response is reduced. This 
value of frequency drop speed can activate the protection operation responding to df/dt. 
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Figure 9. Behavior of the frequency under faults with the presence of the arc 
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Figure 10. Frequency behavior in the event of a two-phase short circuit to ground in 75% of the line that 
interconnects nodes 9 and 10 for different penetration level of photovoltaic solar energy 


4. CONCLUSION 


The studies on frequency instability consider analysis of cases where large load generation 
imbalances occur, either due to a loss of generation or the entry of a load block. The occurrence of short 
circuits usually causes an acceleration of the synchronous units and therefore an increase in the frequency, 
this occurs in systems with large transmission networks, however in systems where networks with an R/X 
values greater than | predominates. In these cases, a short circuit can cause a sudden frequency drop 

There are several factors that condition the response of generators to faults, the type of fault, the 
place where the fault occurs and whether it is with the presence of an arc or not, they affect the module and 
voltage angle and these in the electrical power delivered. The renewable energy sources penetration in the 
composition of the generation further accentuates this phenomenon of the frequency drop as it decreases the 
kinetic energy stored in the system, given the proportion of synchronous units in the total generation and 
therefore before it. Faults are more pronounced and frequency drops at higher speeds. The results showed 
that in weak networks that present units that do not contribute inertia to the network, short-circuit analyzes 
should be included to consider actions in order to prevent the occurrence of frequency instability. 
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